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1. Introduction 



In a recent paper [p we proposed a chiral non-relativistic EFT which included two dibaryon 
fields [Q] as a fundamental degrees of freedom. This EFT, which will be simply referred as NNEFT 
in this paper, is renormalizable and has simple counting rules when dimensional regularization 
(DR) and minimal subtraction (MS) scheme are used. The nucleon-nucleon scattering amplitudes 
in the l So and 3 Si channels were calculated at NLO and a good description of data achieved in 
the — 50MeV energy range. We report here on a calculation at N 2 LO, which was carried out 
in [|3j], in order to see if the good description of data persists and check the convergence of the 
EFT. This is in fact mandatory in view of the fact that the so called KSW approach 0, |||] also 
produced a good description of data at NLO, but turned out to have a bad convergence in the 
3 Si at N 2 LO ^ We will restrict ourselves to an energy range E such that E <^m n , the pion 
mass, and p = ^Jm^E ~ m % . Pion fields can then be integrated out leading to an EFT, which was 
already described in [jlj], which we will call potential NNEFT (pNNEFT). For p ~ m K this EFT 
will already be suitable to carry out the calculations of the amplitudes. For p < 7^ however it will 
be convenient to integrated out nucleon fields with p ~ m K and use the so called pionless NNEFT 
(^NNEFT). All matching calculations will be done expanding the low energy or momentum scales 
in the integrals and using dimensional regularization to regulate any possible UV divergence and 
minimal subtraction scheme. Local field redefinitions which respect the counting will be used to 
get ride off redundant operators, rather than using the on-shell condition. 

2. The nucleon-nucleon chiral effective theory with dibaryon fields 

Our starting point is the effective theory for the Nb=2 sector of QCD for energies much smaller 
than A x (about 2m^, being the nucleon mass) recently proposed in The distinct feature of 
this EFT is that in addition to the usual degrees of freedom for a NNEFT theory, namely nucleons 
and pions, two dibaryon fields, an isovector (D a s ) with quantum numbers x So and an isoscalar (D v ) 
with quantum numbers 3 Si are also included. Since ~ A^, a non-relativistic formulation of 
the nucleon fields is convenient [§]. Chiral symmetry, and its breaking due to the quark masses in 
QCD, constrain the possible interactions of the nucleons and dibaryon fields with the pions. The 
Nb = and Nb = 1 sectors are the usual ones. We will need only the LO lagrangian in the Nb = 
sector and up to the NLO lagrangian in the Nb = 1 sector. 

The Nb = 2 sector consist of terms with (local) two nucleon interactions, dibaryons and 
dibaryon-nucleon interactions. The terms with two nucleon interactions can be removed by lo- 
cal field redefinitions ^ and will not be further considered. The LO terms with dibaryon fields 
and no nucleons in the rest frame of the dibaryons read, 



D](-id Q + S' m )jD^ +Dl(-idv + S' m )jD v , (2.1) 

where D s = D a s T a and S' m ., i = s, v are the dibaryon residual masses, which must be much smaller 
than A x , otherwise the dibaryon should have been integrated out as the remaining resonances 
have. The covariant derivative for the scalar (isovector) dibaryon field is defined as d$D s = doD s + 
\ [[w, dou],D s ]. The NLO pion-dibaryon lagrangian 
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J^ (;? 2) =s x Tr[D s (uJ£^u + u*J£u*)D]} + s 2 Tr[D] (uJt^u + u*JKu*)D s ] + 

+ S3Tr[DlD s uouo] + s^Tr[D]D s UiUi] + s 5 Tr[DluoD s uo] + s 6 ^[dJw;D,w;] + 
+ viDj • b v Tr[u^Jtu^ + w^w] + v 2 D\ • D v Tr[u u ] + v 3 5J •D v rr[w,-w,-] + 
+ v 4 (D^D* + D^' 1 " ) 7> [ Mi M</ ] , 

Si • ~ i = 1, ..,6 and V/, j = 1, ..,4 are low energy constants (LEC). The dibaryon-nucleon 

interactions will also be needed at NLO. They read 



(2.3) 



^ =^ N ^° 2 ^ T2N *) D ^ + f ^{N' T a 2 T 2 T a N)Dl a + 



+ ^(Nh 2 da 2 N*)-D v + ^(N T z 2 (J 2 dN)-Dl 
V2 v 2 



v2 v2 ' v2 



+ -^(7V t t 2 cj 2 5V 2 7V) Dj + -^(V i NU 2 G i G 2 V j N*)Di + -JL(V i N T T 2 o 2 o i V j N)D J v \ 

(2.4) 

with A^,A V - A x 1/2 ,B s ,B v ,B f v - 1/A % . 

3. Dibaryon propagator and counting 

The tree level dibaryon propagator expression ij (-E + S' m . — ir\) gets an important contribu- 
tion to the self-energy due to the interaction with the nucleons as discussed in [jl|] 

i 

72 i = s,v, O- 1 ) 

which is always parametrically larger than the energy E. The size of the residual mass can be 



extracted computing the LO amplitude using the propagator ( p.lj ) and matching the result to the 
effective range expansion, 



^^-T- ' = 5 < v - (3 - 2) 

where at are the scattering lengths of the l So and 3 S\ channels respectively. As a consequence the 
full propagator can be expanded. Moreover equation ( [Ol ) implies that the dibaryon field should 
not be integrated out unless p <C S' m ., instead of E <C 8^. as the tree level expression suggests. Since 

2 

8' m < it should also be kept as an explicit degree of freedom in the so called ^NNEFT. The LO 
expression for the dibaryon field propagator becomes for p (pNNEFT), 

7T~ — i = s ,v, ( 3 - 3 ) 
Afm N p 
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Figure 1: (a) and (b) are order 0(m 3 J 2 /A%) contributions to the dibaryon residual mass, (a) These three diagrams sum up zero by 
Wigner symmetry, (b) Wigner symmetry is violated by insertions of i(—E + <5 m/ ). Naively we would expect these diagrams to be of 
higher order, 0(m 5 J 2 / A 5 J 2 ), but the energy term is enhanced by the radiation pion up to (9(m|/A|). Hence the cross in this diagram 
stands only for an insertion of the energy, (c) Order 0(m|/A|) contributions to the dibaryon residual mass, only the diagrams with the 
potential pion inside the radiation pion loop contribute, (d) N 2 LO contribution to the dibaryon-nucleon vertex. 



and for p < ^ (pNNEFT), 



X 



72 i = s,v. (3.4) 

u rm i 1 % 

The expanded terms will be taken into account through an effective vertex, which we will denote 
with a cross in the dibaryon propagator. Higher order terms in this expansion will be equivalent 
to multiple insertions of this vertex. Furthermore for p » 8 r m . the LO Lagrangian becomes both 
scale and SU (4) (spin-flavor Wigner symmetric) invariant, if the interactions with pions are ne- 
glected [[lOp . Note also that the pNNEFT propagator is suppressed by a m Jl /A x factor respect to 
the pNNEFT propagator. 

Except for the above mentioned contributions to the self-energy of the dibaryon fields, which 
become LO, the calculation can be organized perturbatively in powers of l/A x . Hence one expects 
that any UV divergence arising in higher order calculations will be absorbed in a low energy con- 
stant of a higher dimensional operator built out of nucleon, dibaryon and pion fields (note that the 
linear divergence in the self-energy can be absorbed in 5^.). 



4. Matching to pNNEFT 

For energies E ~ m 2 n /A x <C m % , the pion fields can be integrated out. This integration produces 
nucleon-nucleon potentials and redefinitions of low energy constants. We will follow the strategy 



of [|11[], which was inspired in the formalism developed in [|12|]. 

The dibaryon residual masses get contributions from ( ^ ) and higher loop diagrams involving 
radiation pions, like the ones in fig.[l]b and fig.[j]c, 



4„^;, + 4 ra ,, I1+S2)+ 44(i) 2 (^) 3 + ( 2/i;8;i , 



Note that because of 5^. < ^ the quark mass dependence of 8 m is a leading order effect. 

The dibaryon-nucleon vertices may in principle get 0(m|/A^) from a pion loop, but they turn 
out to vanish, except for those which reduce to iterations of the OPE potential which will already 
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Figure 2: Diagrams contributing to the and 3 S\ partial wave amplitudes. 



be included in the calculations in pNNEFT and must not be considered in the matching. There is 
however a two loop contribution of this order involving radiation pions from the diagram in flgjljd. 

a - a (^M} <4 - 2) 

In the Nb = 1 sector the nucleon mass is redefined by a factor 8m N , however it can be reshuffled 
into 8 m by local field redefinitions. The dibaryon-nucleon interactions remain the same as in (]2J|), 
except for the values of the A t which get modified. 

Finally, in the Nb = 2 sector, two nucleon non-local interactions (potentials) due to the one 
pion exchange are introduced. This is the well known one pion exchange (OPE) potential. 



■Znn =\ j d 3 YN^G a TPN{x l )V a p, pa {x l -x 2 )N^G^ a N{x 2 ), (4.3) 



with, 



gi_ f d 3 q qqqp 
~2flJ (2n)W+ml° 



(4.4) 



5. Matching to pNNEFT 

2 

For p < ^ the calculation must be organized in a different way. This is very much facilitated 
if we integrate out nucleon three momenta of the order of m % first, which leads to the so called 
pionless nucleon-nucleon EFT. The Lagrangian of the Nb = 1 sector of this theory remains the 
same as in pNNEFT. For the Nb = 2 sector the only formal difference from pNNEFT is that the 
non-local potentials ( fO| ) become local and can be organized in powers of p 2 jrrv\. Diagrams in 
Fig.Q containing one (or two) potential pion inside a nucleon bubble will contribute to the dibaryon 
time derivative term as well as the dibaryon residual mass. Contributions to the dibaryon time 
derivative can be reabsorbed by field redefinitions of dibaryon fields, while contributions to the 
residual mass simply redefine it. The derivative and non-derivative dibaryon-nucleon vertex get 
contributions from diagrams containing one (or two) potential pion in the dibaryon-nucleon vertex, 
redefining the LEC B t and A t . Analogous diagrams involving three potential pion, not shown in 
this paper, will also analogously contribute to the matching. The OPE potential in ( |4~4| ) becomes 
0(p 2 /m 2 A 2 ) and hence beyond N 3 LO in the p <^m n region. 
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NLO 



N 2 LO 



Figure 3: Diagrams contributing to the mixing ampli- Figure 4: Diagrams contributing to the 3 D X partial 



tude. 



wave amplitude. 





A s (Mev" 1 ) 


$ ms (MeV) 


p < m i/ A x 


0.0243 


-1.51 


p~m n 


NLO 


0.0302 


-2.25 


N z LO 


0.0251 


-10.5 





A v (Mev- 1 ) 


$ mv (MeV) 


B'JA V (MeV^) 


P < m l/ A x 


0.0265 


7.68 




p~m K 


NLO 


0.0359 


13.4 


VIM -6 


N 2 LO 


0.0110 


45.3 





Table 1: Parameters for the channel. 



Table 2: Parameters for the 3 S\ - 3 D\ channel. 



6. The 1 Sq channel 



The l So amplitude is computed from the diagrams in Fig.^ [J3|l . The phase shifts are plotted in 
Fig.0 and Fig.^ and correspond in the 0-5MeV range to the pNNEFT and in the 5-50MeV range 
to pNNEFT. Parameters have been fit to phase shift data, in the 0-5MeV range in the case of the 
pNNEFT and to 5-20MeV range in the case of the pNNEFT. The reason behind the election to 
fit the pNNEFT to the 5-20MeV data range instead of the more natural 5-50MeV data range is 
that the former delivers a more smoother joining between the pNNEFT and the pNNEFT curves. 
Furthermore the convergence behavior is emphasized. Error bands correspond to the size of the 
next order in the phase shift series. 

Both A s and 8 Ms receive N 2 LO corrections when matching from NNEFT to pNNEFT. If the 
whole expressions for A s and 8 Ms were to be used in the N 2 LO amplitude, higher order terms will 
be introduced. Therefore we will differentiate between 8% LO {A N s LO ) and 8% LO (A^ w ) and we will 

2 s s 

plug them in 8 s > N L0 and 8 S ' NL0 respectively. The values of this parameters have been obtained by 
minimizing the sum of the X 2 functions associated to 8 NL0 and 5^ L0 , with the errors given by the 
size of the next order in the phase shift series (i.e (m n /A x ) 2 and (m n /A x ) 3 respectively). 

The pNNEFT amplitude expressions for N 2 LO (i.e. NLO in the pNNEFT counting) and N 3 LO 
(i.e. N 2 LO in the pNNEFT) are formally identical, however theoretical errors are smaller in the 
later. The only diagrams involved are the tree level one (LO) and the tree level with an insertion of 
a cross (NLO and N 2 LO). The pNNEFT phase shift have been fitted independently of pNNEFT. 
Results for the channel parameters are summarized in Table ||. 



7. The 3 Si- 3 £>i channel 

The 3 S\ amplitude is computed from the diagrams in Fig.^ [Q], the mixing amplitude from the 
diagrams in Fig.^[ and the 3 D\ amplitude from the diagrams in Fig.|] The 3 S\ and 3 D\ phase shifts 
are plotted in Fig]7[ Fig.|| and Fig.[T(] respectively. The mixing angle is plotted in Figj9[ In this 
section we analyze 3 S\- 3 D\ channel. We compare the 3 S\ and 3 D\ phase shifts to data as well as 



6 



Nucleon-Nucleon effective field theory with dibaryon fields 



Jaume Tarrus 





Figure 5: Plot of the NLO expression for the % phase shift. The 
blue line shows the Nijmegen data for the 1 Sq phase shift, while the 
red line corresponds to the fit of our expression. 



Figure 6: Plot of the n 2 lo expression for the phase shift. 
As in the previous figure the blue line shows the Nijmegen data for 
the 1 Sq phase shift. 



A'(deg) 




E(MeV) 



Figure 7: Plot of the NLO expression for the 3 Si phase shift. The 
blue line shows the Nijmegen data for the 3 Si phase shift, while the 
red line corresponds to the fit of our expression. 



A'(deg) 




r E(MeV) 
10 20 30 40 50 

Figure 8: Plot of the n 2 lo expression for the 3 S\ phase shift. 
The blue curve is the Nijmegen data for the 3 Si phase shift, while 
red line corresponds to our expression. 



*(deg) 




Figure 9: Plot of the mixing angle. The blue line shows the Ni- 
jmegen data, the green and red lines the NLO and N 2 LO expression 
respectively. 



<5(deg) 




E(MeV) 



Figure 10: Plot of the 3 Di phase shift. The blue line shows the 
Nijmegen data, the green and red lines the NLO and N 2 LO expres- 
sion respectively. 



the mixing angle. The fitting procedure is analogous to the one used for the 1 Sq> channel. However 
in this channel we have minimized the sum of the % 2 functions associated to S V,NL0 and £ n2lo 
obtaining the values ofA^ LO ,8^° and B' v /A v . If we add the % 2 function associated to 8 v ^ 2lo the 
minimization would not converge. Hence A^ lL0 and 8jf L0 have been obtained in a separated fit 
plugging in the values for the NLO parameters from the first fit. e NL0 , § 3d ^ nlo and, S 3DuN2lo do 
not contain free parameters. The 5 Sl plotted correspond in the 0-5MeV range to the pNNEFT and 
in the 5-50MeV range to pNNEFT, the 8 3 ° l and £ plotted correspond to pNNEFT in the whole 
range. Results for the 3 S\- 3 D\ channel parameters are summarized in Table ||. 

8. Conclusions 

We have calculated the nucleon-nucleon scattering amplitudes for energies smaller than the pion 
mass in the and the 3 S\- 3 D\ channels at N 2 LO in a chiral effective field theory which contains 
dibaryon fields as fundamental degrees of freedom. The large scattering lengths in the x So and 
the 3 5i channels force the dibaryon residual masses to be much smaller than the pion mass. We 
organize the calculation in a sequence of effective theories, which are obtained by sequentially in- 
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tegrating out higher energy and momentum scales. We first integrate out energy scales of the order 
of the pion mass. This leads to an effective theory with dibaryon and nucleon fields, pNNEFT. The 
latter interact through potentials. For three momenta of the order of the pion mass, the scattering 
amplitudes are calculated in this effective theory. For three momenta much smaller than the pion 
mass, it is convenient to further integrate out three momenta of the order of pion mass, which leads 
to the so called pionless NNEFT, and carry out the calculations in the latter. Splitting the calcula- 
tion in this way we can take advantage of the modern techniques of the threshold expansions and 
dimensional regularization so that all integrals only depend on a single scale. There is no need to 
introduce a PDS scheme. The technical complexity of the N 2 LO calculation is similar to the one 
in the KSW scheme, but our final expressions are simpler [||]. The numerical results for the phase 
shifts and mixing angle are also similar to the KSW ones. Hence a good description of the x So 
channel is obtained up to center of mass energies of about 35MeV., but for the 3 5i- 3 Di channel our 
results fail to describe data much before: for the 3 5i phase shift and the mixing angle comparison 
with data becomes bad beyond ISMeV., and for the 3 D\ phase shift it is never good. Particularly 
worrying is the fact that for the 3 S\ and the 3 D\ phase shift the N 2 LO calculation compares worse 
to data than the NLO one. The reasons of this failure can be traced back to the iteration of the OPE 
potential, the first N 2 LO diagram in FigQ which gives a very large contribution. 
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